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Abstract. The existence of a neutrino magnetic moment implies contributions to the neutrino mass 
' via radiative corrections. We derive model-independent "naturalness" upper bounds on the magnetic 

\ moments of Dirac neutrinos, generated by physics above the electroweak scale. The neutrino 

mass receives a contribution from higher order operators, which are renormalized by operators 
responsible for the neutrino magnetic moment. This contribution can be calculated in a model 
independent way. In the absence of fine-tuning, we find that current neutrino mass limits imply 
that |jUv| < 10~ 14 Bohr magnetons. This bound is several orders of magnitude stronger than those 
obtained from solar and reactor neutrino data and astrophysical observations. 
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Current advances in uncovering the pattern of neutrino mass and mixing, lead natu- 
rally to questions about more exotic neutrino properties, such as the magnetic moment, 
/l v . In this paper, we describe how the smallness of the neutrino mass may be used to 
set a strong model-independent limit on the size /i v [1]. Neutrino magnetic moments are 
reviewed in [2], and recent work can be found in [3]. In the Standard Model (SM), 
extended to contain right-handed neutrinos, /i v is non-zero but unobservably small, 
/i v ~ 3 x 10~ 19 [m v /l eV] [4]. Current limits are several orders of magnitude larger, 
so a magnetic moment anywhere near the present limits would certainly be an indication 
of physics beyond the SM. The best laboratory limits arise from neutrino-electron scat- 
tering. The weak and electromagnetic contributions to v — e scattering are comparable 
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V^T~l<r 1( \/— , (1) 

where T is the kinetic energy of the recoiling electron. The present limits derived from 
solar and reactor neutrino experiments are |jU v | ^ 1.5 X 10~ 10 /1b [5] and \fl v \ ^ 0.9 x 
10~ 10 /ig [6] respectively. A more stringent limit can be derived from bounds on energy 
loss in stars, |/i v | < 3 x 10~ 12 jUfi [7]. 

The presence of a non-zero neutrino magnetic moment will necessarily induce a 
correction to the neutrino mass term. (The problem of reconciling a large magnetic 
moment with a small mass has been recognized in the past, and possible methods 
of overcoming this restriction through the use of symmetries are discussed in [8].) 
Assuming that /i v is generated by physics beyond the SM at a scale A, its leading 
contribution to the neutrino mass, 8m v , scales with A as 
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where 8m v is a contribution to m v arising from radiative corrections at one-loop order. 
The A 2 dependence arises from the quadratic divergence of the dimension four neutrino 

mass operator, 6"^ = (L^)vr. Although the precise value of this term cannot be calcu- 
lated in a model-independent way, we can estimate that for A ^ 1 TeV and 8m v < 1 eV, 
we require |/i v | ^5 10 14 /lfi- Given the A 2 dependence, this bound becomes considerably 
more stringent for A well above the electroweak (EW) scale. However, if A is not sig- 
nificantly larger that the EW scale, higher dimension operators are important, and their 
contribution to m v can be calculated in a model independent way. 

We start by constructing the most general operators that could give rise to a mag- 
netic moment operator, Vi0^ v F^ v Vr. Demanding invariance under the SM gauge group 
SU(2)l x U(l)y, we have the following 6D operators 

4 6) = ^L^v^v , 4 6) = glT>~a^v^ v . (3) 

After spontaneous symmetry breaking, both and contribute to the magnetic 
moment. Through renormalization, these operators will also generate a contribution to 
the 6D neutrino mass operator 

®m = ^L$v R (V ) . (4 ) 

The three operators, j , 6$ , 6$ j constitute a closed set under renormalization, so 
that our effective Lagrangian is given by 

^eff = C B (Ai)4 6) +^(M)^ 6) +CM(Ai)^ 6) , (5) 

where the operator coefficients, C ; (/i), depend upon the energy scale \x. The magnetic 
moment and mass are related to the operator coefficients as 
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= -AVl(^pj[C B {v)+C w {v)] , (6) 
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8m v = -Cm(v)^- 2 . (7) 

To connect /i v with m v , we thus need to find the relationship between the coefficients 
Cj(ju) at the weak scale, /l = v. This requires that we solve the renormalization group 
equations (RGE) which relate Q(A) to Cy(v). 

Figures (1,2) display representative examples of the one-loop diagrams which con- 
tribute to the renormalization of the 6D operators. (See Ref. [1] for further details.) 
Solving the RGE, retaining only the leading logarithms, we find that /i v and m v are 
related as 
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where 6 W is the weak mixing angle, 
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/= (l-r)-^rtan 2 %-^(l + r)tan 4 % , (9) 



FIGURE 1. Renormalization of the mass operator, 0^ , due to insertions of (left); self- 

renormalization of (right). 
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FIGURE 2. Self-renormalization of M 
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and r is a ratio of operator coefficients at scale A, r= [Cg(A) — Cw(A)] / [Cg(A) +C^(A)]. 
For A ^ 1 TeV, the bound becomes 

l^gxlO-^xf^l (10) 
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and so for / ~ 1 and m v f$ 1 eV, we find |/i v | ~ 10~ 14 /lfi- In principle, larger values of jl v 
could be obtained, but only by fine-tuning the coefficients Q(A) to arrange cancellations 
in Eq. 9 such that /< 1. We therefore conclude that the natural size of /i v for Dirac 
neutrinos is at least 10 2 times stronger than astrophysical limits and 10 4 times stronger 
than reactor and solar neutrino bounds. The limits that can be placed on transition 
magnetic moments of Majorana neutrinos are substantially weaker than the Dirac case, 
and have recently been calculated in [9]. 
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